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ABSTRACT 


He  have  performed  self-consistent  Linearized 

Augmented  Plane  Have  calculations  of  the  charge  density  of 

* 

metallic  beryllium,  with  and  without  warping  of  the 
interstitial  potential,  using  the  Local  Density 
Approximation  for  the  exchange  and  correlation  energy 
functional.  Overall  agreement  with  the  newest  x-ray  form 
factors  is  good,  and  it  is  argued  that  the  remaining 
discrepancies  have  their  source  in  the  inadequacy  of  the 
Local  Density  Approximation  and,  for  some  reflections,  in 
neglecting  the  non-sphericity  of  the  core  charge  density. 


1.  INTRODUCTION 


Recently  very  careful  measurements  of  the  X-ray 
structure  factors  for  Be  metal  have  become  available  (Hansen 
et  al  1984).  Since  the  X-ray  structure  factors  are  directly 
related  to  the  electron  charge  density,  which  in  principle 
is  exactly  calculable  from  Density  Functional  Theory  (DFT) , 
these  results  can  be  used  as  a  stringent  test  of  present  day 
practice  of  DFT  i.e.  the  Local  Density  Approximation  (LDA) 
for  the  exchange  and  correlation  energy  f uncional ( f or  recent 
reviews  see  Theory  of  Inhomogenuous  Electron  Gas, 
ed . S . Lundqvist  and  N.H. March,  Plenum  Press,  New  York,  1983). 
In  fact  the  X-ray  structure  factors  provide  a  more  severe 
test  of  LDA  than  its  very  successful  application  to  the 
calculation  of  such  physical  quantities  as  cohesive 
energies,  compressibilities,  lattice  constants  etc.,  which 
are  obtained  from  total  energy  differences  where  substantial 
cancellation  of  errors  can  occur  (see  Morruzzi  et  al.  1973). 
No  such  cancellation  of  errors  takes  place  in  the 

I 

calculation  of  the  structure  factors. 

It  is  important  to  be  aware  that  in  the  practical 
implementation  of  DFT  there  are  two  distinct  sources  for  the 
discrepances  between  theory  and  experiment.  One,  which  may 
be  termed  as  fundamental,  is  due  to  the  fact  that  the  exact 
exchange-correlation  energy  functional  is  unknown .  As  in 


most  current  applications,  the  Local  Density  Approximation 
(LDA)  (see  eg.  Williams  and  von  Barth,  1983)  is  used  in  the 
present  work  ( specif ically  the  standard  form  for  exchange 
with  the  Vosko  et  al.  ( 1980 ), hereafter  referred  to  as  VWN, 
form  for  the  correlation  part). The  second  comes  from  the 
computational  aspect  of  the  problem,  namely  the  numerical 
solution  of  the  self-consistent  Kohn-Sharn  equations  for  the 
Belectrons  in  a  crystal,  for  which  a  number  of 
approximations  must  be  made.  Although  in  principle  these  are 
controllable ,  increasing  the  accuracy  of  the  results  puts 
rapidly  growing  demands  on  the  necessary  computer  resources. 
In  the  present  case  we  only  require  that  the  degree  of  self- 
consistency  be  such  that  the  errors  in  the  calculated  form 
factors  be  less  than  the  experimental  ones. 

To  investigate  the  effects  of  the  metallic 
environment  on  the  charge  redistribution  in  Be,  we  have 
performed  Linearized  Augmented  Plane  Wave  (LAPW) 
calculations  for  the  valence  electrons , with  and  without 
warping  of  the  interstitial  potential.  The  Is  electrons  were 
treated  as  atomic-like,  since  they  had  totally  negligible 
amount  of  charge  in  the  interstitial  region  (less  than  0.001 
of  an  electron).  However,  they  were  not  frozen,  i.e.  their 
wave  functions  were  recalculated  at  each  iteration  step. 


We  have  used  experimental  lattice  constants:  a  = 
4.3211  and  c  =  6.7715  a.u.  The  resulting  Fourier  components 
of  the  charge  density  compare  favourably  with  experimental 
data,  as  discussed  below.  However,  a  number  of  discrepancies 
remain.  We  argue,. by  reference  to  atomic  LDA,  Hartree-Fock 
(H-F)  and  Conf iguration  Interaction  (Cl)  charge  densities, 
that  the  discrepancies  for  high  Fourier  components  of  the 
charge  density  are  due  to  inadequacy  of  LDA  to  describe 
electron-  electron  interaction  in  the  core  region,  where 
electron  density  gradients  are  large.  We  also  expect  that 
an  appreciable  discrepancy  for  the  second  Fourier  component 
is  due  to  the  fact  that  we  have  not  taken  into  account 
properly  the  distortion  of  the  Is  core  caused  by  the 
metallic  environment. 

2.  OVERLAPPING  ATOMIC  CHARGE  DENSITIES 

As  a  first  step  we  calculate  the  form  factors  for  the 
density  distribution  corresponding  to  overlapping  Be  atomic 
electron  densities  obtained  using  different  treatments  of 
many-body  effects.  The  intention  is  to  calibrate  the  LDA 
with  respect  to  more  exact  treatments  (i.e.  Cl,  which  we 
consider  virtually  exact,  and  H-F),  to  have  a  reference 
system  from  which  to  measure  metallization  effects  and 


finally  to  have  an  estimate  of  the  relative  importance  of 
core  versus  valence  electron  contributions  to  the  form 
factors . 

The  results  are  summarized  in  Table  I.  The  third 
column  gives  form  factors  obtained  using  atomic  charge 
densities  from  a  Cl  calculation  using  55  configurations 
(Boyd  1979),  the  fourth  column  presents  results  from  LDA, 
and  the  fifth  -  H-F  calculation  using  Clementi  and  Roetti 
(1974)  wave  functions.  For  H-F  results  we  also  show 
separately  the  contribution  from  the  2s  (valence)  electrons 
Finally  column  two  shows  the  experimental  results  of  Hansen 
et  al.(1984),  corrected  for  effects  of  atomic  vibrations 
(Larsen  et  al  1980),  together  with  estimated  errors. 

From  a  comparison  of  the  three  theoretical  columns  it 
can  be  seen  that  the  H-F  values  are  much  closer  to  the 
"exact"  Cl  values  than  the  LDA  ones.  This  is  in  agreement 
with  previous  observations  that  the  local  treatment  of 
exchange  is  a  more  important  source  of  error  in  the  charge 
density  than  the  neglect  of  correlation  (Vosko  and  Wilk 
1983).  Moreover  ,  the  difference  between  H-F  and  Cl  is  of 
the  order  or  less  than  the  experimental  uncertainty, 
suggesting  that  if  one  could  perform  a  calculation  for  the 
metal  equivalent  to  H-F  (e.g.  essentially  exact  treatment  o 
exchange  within  DFT)  the  error  would  be  much  smaller  than 


the  experimental  error  and  the  main  source  of  error  would  be 
in  the  accuracy  of  solving  the  Kohn-Sham  equations.. 

Comparing  Cl  and  the  experimental  results  for  the 
individual  wave  vectors  one  sees  that  in  most  cases  they 
differ  by  more  than  the  uncertainity  in  the  measurements 
indicatincr  the  necessity  of  including  the  effects  of  the 
crystal  environment.  On  the  other  hand  the  magnitude  of  the 
differences  generally  decreases  with  increasing  lenqrth  of 
wave  vector,  as  expected,  since  crystal  field  effects  are 
smaller  in  the  vicinity  of  the  nucleus,  which  becomes  more 
important  as  the  order  of  reflection  increases.  However, 
since  the  LDA  is  used  in  the  metallic  calculation,  it  is 
also  important  to  note  that  the  magnitude  of  the  percentage 
difference  between  the  atomic  Cl  and  LDA  calculations 
increases  from  ~  0.5%  to  ~  1%  as  the  order  of  reflection 
increases.  Thus  we  expect  the  difference  in  the  LDA  metallic 
calculation  and  experiment  to  increase  with  the  magnitude  of 
the  wave  vector.  This  is  in  fact  what  occurs  as  can  be  seen 
from  Table  II. 

The  other  fact  we  would  like  to  point  out  is  that  the 
contribution  of  the  2s  states  is  only  1-2%  of  the  total. 

Thus  even  thouah  we  expect  the  effect  of  the  metallic 
environment  on  the  core  states  to  be  small,  it  can  have  an 
appreciable  influence  on  the  chancre  of  form  factors  on 


goincr  from  atom  to  metal,  since  the  contribution  from  the 
core  is  the  dominant  part  of  the  total  form  factor. 


3.  LAPW  FORM  FACTORS. 


3.1  Calculational  details 


To  study  effects  of  metallization  on  the  electronic 
charge  distribution  we  have  performed  self-consistent, 
warped  muffin-tin,  LAPW  band  structure  calculations  (see 
e.g.  Andersen  1975,  Koelling  and  Arbman  1975),  with  the  VWN 
form  of  LDA.  Basis  functions  inside  the  muffin-tin  included 
anpular  momentum  components  up  to  Jl  =  8  and  in  the 
interstitial  region  we  found  it  necessary  to  go  to  over  50 
plane  waves,  i.e.  to  take  k.  'Rwrr,  -  0.4.  Increasing 
^max’^MT  7  •  ^  did  not  modify  the  results  in  a  significant 
way.  The  energy  parameter  was  chosen  in  the  middle  of  the 
band,  E=0.4  Ry.  Brillouin  zone  integrations  were  performed 
using  the  special  points  method  (Fehlner  and  Vosko  1977), 
and  we  used  as  a  standard  35  points  in  the  irreducible  part 
of  the  B.Z.  Increasing  this  number  to  57  had  only  a  minor 
effect  on  the  form  factors,  rionspher icity  of  the  valence 
charge  density  inside  the  muffin-tin  was  taken  into  account 
only  approximately  by  analytically  continuing  the  non- 
spherical  part  of  the  interstitial  charge  density  into  the 


pseudopotential  method.  The  MT  calculation  was  performed  to 
investigate  the  importance  of  "warping"  of  the  interstitial 
potential  on  the  charge  density.  For  the  low  order  beams 
(with  the  exception  of  002)  the  warped  calculation  agrees 
better  with  experiment  than  the  atomic  calculations.  Warping 
is  important  for  obtaining  agreement  for  the  lowest  order 
reflection  but  affects  the  other  reflections  only  slightly. 
The  disagreement  in  the  002  form  factor  is  probably  due  to 
our  neglect  of  the  non-spherici ty  of  the  core.  Since  the  c/a 
ratio  (1.57)  for  Be  metal  is  appreciably  different  from  that 
for  an  ideal  hep  crystal  (1.63)  it  is  likely  that  the  core 
charge  would  be  more  non-spherical  than  in  a  close-packed 
structure.  The  002  component  is  sensitive  only  to 
distortions  in  the  z  direction. 

Table  TII  compares  the  form  factors  of  the  core  charge 
density  from  our  self-consistent  calculations  with  those  for 
a  free  Be  atom  using  the  VWN  form  of  LDA.  The  difference 
between  these  form  factors  is  shown  in  the  fourth  column  and 
the  difference  divided  by  the  corresponding  geometric 
structure  factor  in  the  fifth  column.  The  entries  in  the 
fifth  column  are  all  positive  and  increase  with  the 
magnitude  of  sinO/X,  indicating  an  expansion  of  the  cores  in 
the  metallic  environment  as  expected. 


The  values  of  the  last  three  form  factors  obtained  from 


the  warped  calculation  are  all  appreciably  smaller  than  the 
experimental  values.  This  is  probably  due  to  the  inadequacy 
of  LDA  when  the  gradients  of  the  charge  density  are  large. 

_  The  LDA  is  known  to  give  much  better  results  when  used 
to  calculate  differences  of  various  quantities  than  their 
absolute  values.  We  can  crudely  include  both  the  effects  of 
metallization  and  non-local  effects  in  the  exchange 
correlation  potential  by  recasting  our  results  in  terms  of  a 
difference  in  the  two  LDA  calculations  plus  the  Cl  form 
factor  (using  H-F  atomic  form  factor,  if  Cl  was  not 
available,  would  also  be  adequate): 


fD  =  fLDA(met)  ~  fLDA(at)  +  fCI(at) 


(1) 


where  fLD^(met)  is  our  calculated  form  factor  for  the  metal 
(column  2  of  Table  II),  is  atomic  form  factor 

(column  4  of  Table  I)  -  both  calculated  within  LDA  and 
f^-j-(at)  is  the  Cl  atomic  form  factor  (column  3  of  Table  I). 
We  see  that  in  creneral  we  qet  an  improvement ,  with  the 
exception  of  the  002  reflection,  which  becomes  worse.  This 
again  suggests  that  it  is  distortion  of  the  cores  due  to  the 
metallic  environment,  which  is  a  source  of  the  discrepancy 
in  this  case. 


4.  CONCLUSIONS 


In  conclusion  we  stats  that  metallization  effects  are 
important  for  obtaining  agreement  of  theoretical  short  wave 
vector  Fourier  components  of  charge  density  in  metallic  Be 
with  recent  accurate  experimental  data.  In  particular,  for 
the  shortest  component  inclusion  of  warping  of  the 
interstitial  potential  is  essential.  For  higher  q  components 
the  LDA  seems  to  become  less  adequate,  which  is  not 
surprising  since  such  Fourier  components  are  sensitive  to 
the  charge  distribution  in  the  core  region  where  gradients 
of  the  charge  density  are  large  and  hence  corrections  to  LDA 
are  expected  to  be  substantial.  We  also  note  that  the 
contribution  from  valence  electrons  constitutes  less  than 
4%  of  the  total  value  of  the  form  factors  (and  this 
fraction  goes  down  with  increasing  q)  and  hence  even  a  small 
distortion  of  the  core  wave  functions  may  contribute 
significantly  to  the  form  factors.  We  expect  that  this 
effect,  not  taken  into  account  in  our  calculations,  is 
responsible  for  marked  disagreement  between  our  calculated 
and  experimental  values  for  the  002  form  factor. 
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TABLE  CAPTIONS 


Table  I 

Form  factors  of  overlapping  atomic  charge  densities  for  Be, 
calculated  within  different  approximations: 

Cl  -  Configuration  interaction  with  55  configurations  (Boyd 
1979;  ) 

LDA  -  Local  Density  Approximation  in  VWN  form  (Vosko  et  al 
1980  )  ; 

H-F  -  Hartree-Fock.  using  Clementi  and  Roetti  (  1974)  wave 
functions ; 

H-F  2s  (column  6)  is  the  contribution  from  H-F  ws  wave 
function.  Experimental  data  are  from  Hansen  et  al.(1984) 

Table  II 

Form  factors  of  charge  density  for  metallic  beryllium. 
Columns  2  and  3  show  our  LAPN  results  with  (warped)  and 
without  (MT)  warping  of  the  interstitial  potential.  For 
comparison  we  also  show  the  results  of  Chou  et  al.(1983)  in 


column  4.  Experimental  data  are  from  Hansen  et  al.<1984). 
Column  6  gives  the  values  of  the  quantity  fQ  defined  in  the 
eq.  1 . 

Table  III 

Structure  factors  of  the  core  charge  density  for  Be.  The 
second  column  gives  values  obtained  in  our  "warped"  LAPW 
calculation,  the  third  column  refers  to  the  LDA  atomic 
calculation,  the  fourth  column  is  the  difference  of  the 
previous  two  and  in  the  last  one  we  give  this  difference 
divided  by  the  geometric  structure  factor  (G5F) . 
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1.817 

1.818 

0.020 

002 

3.432(7) 

3.491 
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3.482 

-0.034 
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2.938(7) 
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110 

2.909(18) 

2.921 
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-0.036 
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2.391 

-0.006 

200 

1.368(7) 

1.351 

1.337 

1.353 

0.002 

004 

2.506(10) 

2.527 

2.497 

2.531 

0.024 

TABLE 

II 

Present 

calc . 

:1  warped 

MT  . 

Chou  et  al 

Exp . 

fD 

)0  1.886 

1.846 

1.8845 

1.889(6) 

1.895 

)2  3.468 

3.465 

3.4418 

3.432(7) 

3.486 

11  2.929 

2.944 

2.9255 

2.938(7) 

2.944 

>2  1.550 

1.558 

1.5365 

1.550(11) 

1.561 

.0  2.889 
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13  2.348 

2.355 

2.3362 
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2.372 

'0  1.322 

1.329 

1.316 
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1.336 
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2.478 

2.4798 

2.506(10) 
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